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Abstract 

The meson-photon transition form factors 77* — > P {P stands for tt, n and rf) provide strong 
constraints on the distribution amplitudes of the pseudoscalar mesons. In this paper, these tran- 
sition form factors are calculated under the light-cone perturbative QCD approach, in which both 



the valence and non-valence quarks' contributions have been taken into consideration. To be con- 
sistent, an unified wavefunction model is adopted to analyze these form factors. It is shown that 



> 

IT) 

I with proper charm component fZ ~ —30 MeV and a moderate DA with B ~ 0.30, the experi- 

m ■ 

mental data on Q 2 F ril (Q 2 ) and Q 2 F„ij(Q 2 ) in whole Q 2 region can be explained simultaneously. 



Further more, a detailed discussion on the form factors' uncertainties caused by the constituent 
quark masses m q and m s , the parameter B, the mixing angle <p and f^, are presented. It is found 
that by adjusting these parameters within their reasonable regions, one can improve the form factor 
to a certain degree but can not solve the puzzle for Q 2 F 7T ~ f (Q 2 ), especially to explain the behavior 
of 7r — 7 form factor within the whole Q 2 region consistently. We hope further experimental data 
on these form factors in the large Q 2 region can clarify the present situation. 
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I. INTRODUCTION 



The distribution amplitude (DA) is the key factor for exclusive processes. Usually, the 
DAs of the light pseudoscalar mesons can be expanded in Gegenbauer polynomials, and their 
corresponding Gegenbauer moments have been studied by various groups, cf. Refs.[l-9]. 
However, there is no definite conclusion on whether the shape of the DA is in an asymptotic 
form [lOj] or in a more broad form 11]. 

The light pseudoscalar meson-photon transition form factor Fpj(Q 2 ) that describes the 
effect of the strong interaction on 77* — > P transition (P stands for 71", 1] and 7/) provides 
a good platform for studying the leading-twist DA, since it contains only one bound state 
and the power suppressed light meson's higher helicity and higher twist structures usually 
give negligible contributions. One can extract useful information on the leading-twist DA 
by comparing the theoretical with the measured data on Fp^(Q 2 ). 

Based on CELLO, CLEO and the BABAR data on 77* — > 11 [rj-yj], many people have 
discussed the properties of pion DA 1514281] . A CZ-like DA or even a flat DA [15J can explain 
the large Q 2 behavior shown by the BABAR data [14], however the theoretical estimation on 
the form factor with a CZ-like or flat DA shall always be lower than the experimental results 
in small Q 2 region 22]. While by taking the non- valence quark parts into consideration, it 
is found that by setting the second pion moment a^/i 2 ,) around 0.35, one can explain the 
behavior in small Q 2 region well, however a somewhat large discrepancy emerges in high 
Q 2 region if the BABAR data is confirmed 28J]. A reasonable theoretical estimation on the 
form factor should explain the measured form factor's behaviors in both the lower and higher 
Q 2 -regions consistently. And one should find a way to compare the experiment results on 
the form factor to determine which DA shape is more suitable. 

More over, the experimental data on 77* — > 77 or rj transition 
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35|, especially, 



the new BABAR results within the reg ion of [4, 40] GeV 2 36|, can provide further constraints 
on the pseudoscalar meson's DA 37H43|. Then, we shall have three pseudoscalar meson- 
photon transition form factors Fp 7 (Q 2 ) to constrain the light meson's wavefunction (WF) / 
DA parameters. Because of 77 and 7/ mixing, their condition is somewhat more difficult than 
the pionic case. Two mixing schemes are adopted in the literature, even though their mixing 



parameters can be related through the correlation given by Refs. 



44 



45] , there are differences 



in dealing with certain processes. One mixing scheme is based on the flavor singlet i]i and 
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octet 7/s, under which, one usually introduces two mixing angles 6\ and 6 8 and adopts the 
same DA for r/iand rjs 38|-|43| . In the present paper, we adopt the simpler quark-flavor 
mixing scheme [441 ]. which is based on the quark-flavor basis r\ q and r\ s and only one mixing 
angle <fi is introduced. Since r\ q and 7] s have similar structure as that of pion, it is natural to 
adopt the same WF model for 7r, r\ q and r\ s . For the purpose, we adopt the WF model raised 



by Ref . 



281 ] to do our discussion, since by setting the parameter B properly, we can obtain 



different DA behavior from AS-like to CZ-like naturally, and then to determine which one 
is more suitable for simultaneously explaining the experimental data of these form factors. 

The paper are organized as follows. In Sec. II, we outline our calculation techniques for 
the transition form factors -Fp 7 (<5 2 ), where the mixing scheme for 77 and rj', an uniform WF 
model for the mesons, and the analytic formulas for deriving Fp 7 (Q 2 ) are presented. In 
Sec. Ill, we present the numerical results, and discuss the uncertainty sources. The final 
section is reserved for a summary. 



II. BASIC FORMULAS FOR THE FORM FACTORS F Py (Q 2 ) 

In this section, we present necessary formulas for the transition form factors Fp^(Q 2 ). 
First, we define the physical meson states 77 and 7/ under the quark-flavor basis, then we 
give the uniform WF model for the mentioned pseudoscalar mesons, and finally, we present 
the form factor with the valence quark contribution calculated up to next-to-leading order 
(NLO) together with an estimation of the non-valence quark state's contributions. 



A. r] and 7/ denned under the quark-flavor basis 

The physical meson states rj and 7/ are related to the orthogonal states r\ q and r\ s through 
an orthogonal transformation 

/ \v } \ f\Vq)\ /cos0 -sin0\ 

\\V )/ \\Vs) / \sm0 cos0/ 

where <fi is the mixing angle. Here, we adopt a single mixing angle scheme that attributes 



SUf(3) breaking to the Okubo-Zweig-Iizuka violating contribution 44|. In the quark-flavor 
basis, the two orthogonal states \r) q ) and \r) s ) are defined in a Fock state description, \r] q ) = 
^ \um+Mi and |^ = vj/^Jss^ where ^ (i = q, s) denote the light-cone WFs of the 
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corresponding parton states. Under such scheme, the decay constants in the quark-flavor 
basis simply follow the pattern of state mixing, i.e. 

f ] f ) = U(<j>)dmg[f Vq J Vs ], (2) 

where f Vi = (2a/3) / fe 2 <fJ 2 ^fjt- ^ m (x, k±) and the factorization scale /x ~ 0(1 GeV). 

One can obtain the correlation between f Vq /f Vs and (f) from the two-photon decay of rj 
and if , i.e. 77/77' — )■ 77, which shows [37] 



and 
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(4) 



where c s = y/2/3, c q = 5/3 and a = 1/137. Since the power suppressed higher twist and 
higher helicity components give negligible contributions to the two-photon decay of 77 and 
r]', the correlations ( I3f4"l) shall provide strong constraint on f Vq , f Vs and (j). 



B. wavefunction of the light pseudoscalar meson 



As for the light pseudoscalar meson (P), its light-cone WF can be written as 
* P (x,k ± ) = X AlA2 (2^,k±)^p(^™*,k ± ), 

Ai A2 



(5) 



where i stands for the light quark q or s, Ai and A2 are helicities of the two constituent 
quarks. x AlA2 (x, kj_) stands for the spin-space WF coming from the Wigner-Melosh rotation. 
*&fq(x, rrii, kj_) stands for the spatial WF, which can be factorized as [28] 



typ(x, mi, kj_) = Aipp{x) exp 



lpi 2 x(l - x) 



(6) 



The x-dependent (pp(x) can be expanded in Gegenbauer polynomials, and by keeping its 
first two terms, we obtain 



m h kj.) = A(l + B x Cl /2 {2x 



1)) exp 



+ to- 
3/3 2 x(l -x) 



(7) 
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As for the parameters B and (3, B determines the broadness of the WF in the longitudinal 
direction, while (3 determines the WF's transverse behavior. 

The DA 4>p{x) can be obtained from typ(x,kj_) by integrating over the transverse mo- 
mentum, <pp(x) = J| P 



|k x | <w , 1& (7?) M^), and we obtain 



(f> P (x,i4) 



VSAm/3 



x(l-x)[l + Bx Cl ,2 (2x- 1) 



^Erf 




m 2 + ill 


-Erf 








%f3 2 x(l - x) 







\ S(3 2 x(l - x) 



(8) 



where Erf(x) 

nomials, whose moments a n (/j, 2 



^= Jq e t2 dt. (f) P (x, Hq) can be expanded in conventional Gegenbauer poly- 



Jl dx^pjx^Dcl' 2 \2x-l) 



Numerically, it is found that its 



° J Jl dx6x(l-x)[C 3 J 2 (2x-l)}2 

second Gegenbauer moment 02(^0) * s c l° se to B, i.e. the DA's behavior is dominated by 
B. Moreover, when B ~ 0.00, its DA is asymptotic-like; and when B ~ 0.60, its DA is 
CZ-like. This shows (f)p(x, /jLq) can mimic the DA behavior from asymptotic-like to CZ-like 
naturally. Then, by comparing the estimations for B e [0.00, 0.60] with the experimental 
data on various processes, one may decide which is the right DA behavior possessed by the 
light pseudoscalar mesons. Here, we do not discuss the flat DA, since it is hard to explain 
the meson-photon transition form factor's behavior around Q 2 ~ and shall meet even more 
serious end-point problem than the CZ-like DA at x ~ 0, 1 



C. Pseudoscalar-photon transition form factors 



The pseudoscalar meson-photon transition form factors can be divided into two parts 



Fpi{ Q'2 )=F (V) {Q 



i(NV) 



(Q 5 



(9) 



P 7 IV / ' P*y 

where F P v y \Q 2 ) is the valence quark part, F p n 1 v \Q 2 ) stands for the non- valence quark part 
that is related to the higher Fock state of the pseudoscalar meson. 

Under the light-cone perturbative QCD approach 10| , and by keeping the /^-corrections 
in both the hard-scattering amplitude and the WF, F n7 (Q 2 ) has been calculated up to NLO 
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511 ] . For pseudoscalar meson-photon transition form factors up to NLO, we have [28[ 



f { pV(Q 2 



e P r 1 r x2 Q 2 dx 
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(10) 



where [dx] = dxdx'5(l — x — x'), Cf = 4/3 and k± = |kj_|. And ep relates to the electric 
charge of the constituent quarks, e T = 1/3, e Vq = 5/9 and = a/2/9. 

As for Fp^ V \Q 2 ), we adopt the model suggested by Refs. 47], which is constructed based 
on the form factor's limiting behavior at both Q 2 — > and Q 2 — > oo, i.e., 



7 (NV) 

Pi 



(Q 2 ) 



a: 



(11) 



where a = |Fp 7 (0), k 



pi 



W 2 )I Q 2^0 



(l + Q 2 //€ 2 ) 2 ' 

with the first derivative of F p ^ V \q 2 ) over Q 2 



in the limit Q 2 — > takes the form 



(Q 2 



e P 



8tt 2 



d 



1 /-x 2 Q 2 



<9Q 2 "'O Jo 



Q 2 ^0 



The same phenomenological model for F p N y V \Q 2 ) has also been adopted by Ref.[l6|, 

GeV is introduced to replace the parameter 



where instead a fixed input parameter A ~ 1. 
k. Since the octet-singlet mixing scheme 44, 



45| is adopted by Ref.|l6|, it is reasonable to 



take the same A for n, r]g and t]q. Under the present adopted quark-flavor mixing scheme, 
we shall numerically obtain k ~ 1.1 — 1.2 GeV for pion and r] q , and k ~ 1.5 — 1.6 GeV for 
r/ s , where different k is rightly caused by S , f//(3)-breaking effect. 

Moreover, under the quark-flavor mixing scheme, rj — 7 and 7/ — 7 transition form factors 
are related with F Vq ^(Q 2 ) and F Val (Q 2 ) through the following equations 



F m {Q 2 ) = F %7 (Q 2 )cos0-F, s7 (Q 2 )sin0 



(12) 



and 



i^ 7 (Q 2 ) = F ??97 (Q 2 )sin0 + F r , s7 (Q 2 )cos0. 



(13) 



III. NUMERICAL RESULTS AND DISCUSSIONS 



A. Input parameters 



Two photon decay widths of 77 and rf , and their masses can be found in PDG 52| 



IV^ = 0.510 ± 0.026 KeV, M v = 547.853 ± 0.024 MeV, 
r\^ 77 = 4.28 ± 0.19 KeV, M n > = 957.78 ± 0.06 MeV, 
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TABLE I: Typical WF parameters for m q = 0.30 GeV and m s = 0.45 GeV, where cp = 39.5° is 
adopted. 



B\m 


Ar(GeV) 


^(GeV- 1 ) 


^(GeV- 1 ) 


^(GeV- 1 ) 


0.00 \m q 


0.586 


25.06 


26.81 


/ 


0.30|m g 


0.668 


20.26 


21.67 


/ 


0.60|m g 


0.745 


16.62 


17.78 


/ 


0.00|m s 


0.464 


42.23 


/ 


60.58 


0.30|m s 


0.504 


36.97 


/ 


49.58 


0.60|m s 


0.552 


31.24 


/ 


40.72 



and U = 92.4 ± 0.25 MeV. 

A weighted average of the experimental values shown in Ref. 53|, together with two 
experimental values <p = 38.8° ± 1.2° [sj and = 41.2° ± 1.1° [5^, yields = 39.5° ± 0.5°. 
Then, with the help of Eqs.( )3ll4|) . we obtain 



and 



L = 1.07 ±0.03 



/„.//* = 1.44 ±0.08, 



(14) 



(15) 



53|]. 



which are consistent with the "phenomenological" values, f v ~ f n and f Va ~ 1.36/^ 

As for pion WF, its parameters can be determined by its normalization condition and the 
constraint from tiq — > 77. A detailed determination on pion WF parameters can be found 
in Ref. 28|, where (3 n for specified B and quark mass is determined by 



f 2 
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(16) 



Jo 1 tfe* w (ar, k± = 0) 

Because rj q and rj s have similar behaviors as that of 7T, for clarity, we take (3 q = fi-n\m q and 
0s — (3n\m s - Under the condition of B = 0.00, 0.30 and 0.60, typical parameters for the DAs 
of 7T, rj q and rj s are presented in TAB.(jI]), where the mixing angle <fi is fixed to be 39.5°. 

It is noted that by varying B within the region of ~ [0.00, 0.60], the DAs shall vary from 
asymptotic-like to CZ-like form. To show this point more clearly, we draw (f) v (x) and (f) r)s (x) 
in Fig. ([1]), where B = 0.00, 0.30 and 0.60 respectively. 



FIG. 1: DA model © for cf) Vq (x) (Left) and (f) Vs (x) (Right) with fj, = 1 GeV, where B = 0.00, 0.30 
and 0.60 respectively. 

B. Basic numerical results 

Numerically, one may observe that in the large Q 2 region, the leading valence Fock- 
state contribution dominates the form factor Q 2 Fp 1 (Q 2 ), while the non- valence quark part 
Q 2 FpV (Q 2 ) is power suppressed and is quite small, so it is usually neglected in the literature. 
However, Q 2 Fp^(Q 2 ) can provide sizable contributions in the low and intermediate energy 
regions, so one should take it into consideration to make a more sound estimation in the 
whole energy regions. 

As shown above, the parameter B in the unified WF model ([5]) determines the DA 
behavior of the light pseudoscalar mesons. Then, by comparing with the experimental data 
on the pseudoscalar meson-photon transition form factors, it provides us an opportunity to 
discuss the DA properties in a more consistent way. 

Firstly, we present the meson-photon transition form factors within a wide region of 
B e [0.00, 0.601 so as to show which DA behavior is more suitable to explain the data, 

n 

especially the BABAR data [14J. In doing the numerical calculation, we take all the other 
input parameters to be their center values, i.e. m q = 0.30 GeV, m s = 0.45 GeV and 
(f) = 39.5°. Figs. ([2} [3]) show the pseudoscalar-photon transition form factors Q 2 F nl (Q 2 ), 
Q 2 F m (Q 2 ) and Q 2 F v r 1 (Q 2 ), where B = 0.00, 0.30 and 0.60 respectively. These two figures 
show that with the increment of B, all the three form factors decrease in lower Q 2 region but 
increase in higher Q 2 region. Especially, the CZ-like DA (B = 0.60) leads to the smallest 
value in the lower Q 2 region, while the AS-like one (B = 0.00) leads to biggest one; and 
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FIG. 2: 7r — 7 transition form factor Q 2 F w ~ f (Q 2 ) with varying B. The solid, the dotted and the 
dashed lines are for B = 0.00, 0.30 and 0.60 respectively. The experimental data are taken from 



Refs 



14] 




> 0.4 

O 

o 

LL 0.2 
CM 

o 



• GLEO 
■ BaBar 


data 
data 

t 


** 

* 




. ■ 1 1 ' ** 





Q 2 (GeV 2 ) 



Q 2 (GeV 2 ) 



FIG. 3: r/ - 7 (Left) and 7/ - 7 (Right) transition form factors Q 2 F V1 {Q 2 ) and Q 2 F v/l (Q 2 ). The 



solid, the dotted and the dashed 
data are taken from Refs. 
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33, 



ines are for B = 0.00, 0.30 and 0.60 respectively. The experimental 



36]. 



in the higher Q 2 region, the condition is vice versa. This causes the present puzzle for 
explaining the newly obtained BABAR data on ix — 7 form factor. The BABAR data shows 
that in the range of Q 2 G [4,40] GeV 2 , the pion-photon transition form factor behaves as 

H, Q 2 F nl (Q 2 ) = A (y£y?Y> where A = °- 182 ± °- 002 and f 3 = °- 25 ± °- 02 - A CZ-like 
DA with B ~ 0.60 or even flat DA can explain the data well for higher Q 2 region 1 , however 



1 The form factor with flat DA shall lead to logarithmic growth with Q 2 , i.e. Q 2 F 7r ~ / (Q 2 ) oc In (l + Q 2 /M 2 ) 
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Q 2 (GeV 2 ) 



FIG. 4: 7T — 7 transition form factor Q 2 F 7T ^ f (Q 2 ) with fixed B = 0.3. The red solid line, the blue 



dotted line and the black dashed line are for m 



The experimental data are taken from Refs. 
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0.25 GeV, 0.30 GeV and 0.35 GeV respectively. 



it fails in lower Q 2 region. The AS-like DA with B ~ 0.00 provides a better understanding 
for lower Q 2 region, however due to the fact that Q 2 F ni (Q 2 ) tends to be a constant value 
(2/^), it can not explain the large Q 2 behavior. On the one hand, by increasing B from 
to a larger value, the estimated large Q 2 -behavior of Q 2 F^{Q 2 ) can be improved. On the 
other hand, the deviation of the lower Q 2 -behavior also increases with the increment of B, 
so B should not be too large. Moreover, as shown by Fig. ([3]), rj — 7 form factor Q 2 F V1 (Q 2 ) 
prefers a DA with smaller B, i.e. B < 0.30. For the rj — 7 form factor Q 2 F r f 7 (Q 2 ), all DAs 
shall lead to large Q 2 behavior well above the experimental data. 

Secondly, we study the uncertainties of the transition form factors caused by m q and m s . 
For the purpose, we take m q = 0.30 ± 0.05 GeV and m s = 0.45 ± 0.05 GeV, and fix the 
parameter B = 0.30 and <fi = 39.5°. The 71 — 7 form factor for m q = 0.30 ± 0.05 GeV is 
presented in Fig.fllJ), and rj — 7 and rj — 7 form factors for [m q = 0.25 GeV and m s = 0.40 
GeV], [m q = 0.30 GeV and m s = 0.45 GeV], and [m q = 0.35 GeV and m s = 0.50 GeV] are 
presented in Fig.(JS]). It can be found that the form factors change with the constituent quark 
mass similar to its change with B, i.e. with the increment of constituent quark mass, all 
the three form factors increase in the lower Q 2 region and decrease in the higher Q 2 region. 
Naively, one may expect to obtain a larger high Q 2 behavior by setting a smaller m q and m s . 
Especially, by setting the limiting values of m q = and B = 0, which rightly corresponds to 



22| , which is close to B AB AR data. 
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FIG. 5: 7] - 7 (Left) and ij — 7 (Right) transition form factors Q 2 F m (Q 2 ) and Q 2 F v > y (Q 2 ) with 
fixed B = 0.3 and = 39.5°. The red solid line, the blue dotted line and the black dashed line are 
for [m q = 0.25 GeV and m s = 0.40 GeV], [m q = 0.30 GeV and m s = 0.45 GeV], andjmg ^,0.35 



GeV and m s = 0.50 GeV] respectively. The experimental data are taken from Refs. 



33 



36]. 
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FIG. 6: 7] - 7 (Left) and rj' — 7 (Right) transition form factors Q 2 F m (Q 2 ) and Q 2 F r) >~ / (Q 2 ) with 
fixed B = 0.3 and m q = 0.30 GeV and m s = 0.45 GeV. The experimental data are taken from 



Refs. 



13, 



33, 



36] 



a flat (p n as suggested by Ref. 22| . one can obtain the same logarithmic behavior for large Q 2 



region that is consistent with BABAR data, i.e. Q 2 F 1T1 (Q 2 ) oc \n(Q 2 /2a 2 ) with a = ^-e lE . 
However, it is found that m q can not be too small, e.g. it should be larger than 0.22 GeV, 
otherwise the probability of leading valence quark state \qq) shall be larger that 1 [55]. 

The conditions for 77 — 7 and 7/ — 7 are somewhat more different. Due to i] — 7/ mixing, 
we need to consider these two form factors simultaneously. The curves for Q 2 F rjl {Q 2 ) and 
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Q 2 F v i 1 {Q 2 ) for = 39.5°±0.5° are presented in Fig.®. It is found that Q 2 F m (Q 2 ) decreases 
with the increment of <ft, while Q 2 F V /^(Q 2 ) increases with the increment of 0. By shifting 
4> to a smaller value ~ 38°, one can explain these two form factors within Q 2 < 20 GeV 2 



consistently as has been pointed out by Ref. 37J. However such a shifting of <p can n °t 
explain the newly BABAR data on rj — 7 and rf — 7 for even larger Q 2 > 20 GeV 2 . So 
it is hard to fit the gap between the theoretical estimation and the experimental data in 
the whole Q 2 region by a simple variation of 0. Experimentally, Q 2 F m (Q 2 ) still increases 
with the increment of Q 2 up to a large value even though its ascending trends is slower 
than the growth of Q 2 F WJ (Q 2 ) , while Q 2 F V / 7 (Q 2 ) tends to be a consistent for Q 2 — > 00. So 
some other sources have to be introduced to explain both the 77 — 7 and 77' — 7 form factors 
in the whole energy region consistently As shown by Fig. ([3]), Q 2 F rjl (Q 2 ) can agree with 
the data with B ~ 0.30, so we hope the new sources shall have less effects to Q 2 F V1 (Q 2 ) 
than that of Q 2 F V * 1 (Q 2 ). It has been suggested that a proper intrinsic charm component 



may have some help to explain the abnormally large production of 7/ 53|, l56M58j . In the 
following subsection, we shall make a detailed discussion on the possible contributions from 
the intrinsic charm components. 



C. Possible contributions from the intrinsic charm components to Q 2 F rjl {Q 2 ) and 

Q 2 F V > J (Q 2 ) 



Since the mixing between the cc state with qq-ss basis is quite small |53|], we can set 



(17) 
(18) 



F m (Q 2 ) = F W (Q 2 ) cos <t>-F^{Q 2 ) sin 4, + F^jQ 2 ) 
F V ,,(Q 2 ) = F Vq ,(Q 2 ) sin + F^Q 2 ) cos <p + F^(Q 2 ), 

where F^(Q 2 ) and F^(Q 2 ) corresponds to the contributions from the intrinsic charm 
component in 77 and 77' respectively. Similarly, the WF of the "intrinsic" charm component 
7? c = |cc) can be modeled as 

/ k 2 | + m 2 



% h ,(x, kj.) = A c v/V , (l + B x Cl /2 (2x - r 



cxp 



ip 2 x(l -x)_ 



X (m c ,x,kj 



(19) 



where we adopt (3 C = f}n\m a - The overall factor A^^, is determined by the WF nor- 
malization, in which their corresponding decay constants f° and f5 are related by [531 ] . 



— tan 



— arctan 



Here to calculate F^(Q 2 ) and F^(Q 2 ), the charm quark 
12 
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FIG. 7: 7] - 7 (Left) and ij — 7 (Right) transition form factors Q 2 F m (Q 2 ) and Q 2 F v > y (Q 2 ) with 
varying fZ, where we take i? = 0.30, m q = 0.30 GeV, m s = 0.45 GeV and m c = 1.50 GeV. The 



experimental data are taken from Refs. 



13, 



33 



36]. 



mass effect should be taken into consideration in the hard part of the amplitude, i.e. the 
higher helicity states that are proportional to the quark mass shall provide sizable contri- 
butions. After integrating over the azimuth angle, a direct calculation shows 37] 



where z 



1 dx r°° 



V2 



3\/37T z ./'I •/• -/n 



% /vl (x,k 2 ± ) 



1 



1 - z - y 2 



(z + (l-y)a)(z+(l + y)3) 



k 1 dk 1 



(20) 



and y 



x 2 Q 2 

Taking 5 = 0.30, m q = 0.30 GeV, m s = 0.45 GeV, m c = 1.50 GeV and <p = 39.5°, we 
show how f5 aflects the form factors Q 2 F m (Q 2 ) and Q 2 F V / 1 (Q 2 ). The results are presented 
in Fig.flTJ), where fZ = 0, —5 MeV, —15 MeV and —45 MeV respectively. These two form 
factors are slightly affected by the charm component in low Q 2 region, while in high Q 2 
region, the form factors are quite sensitive to fZ and they can be greatly suppressed by 
possible charm component. One may observe that the experimental data disfavors a larger 
portion of charm component as |/^,| ^ 50 MeV. And for a more larger \fZ\, it shall have 
more obvious effects to Q 2 F rj i^(Q 2 ) than that of Q 2 F vy (Q 2 ), which is what we wanted. 

More over, we present the results for rj — 7 and rj' — 7 transition form factors with 
fixed JZ = —30 MeV in Fig. (JSJ) , where B = 0.0, 0.30 and 0.60 respectively. It shows that 
with proper charm component f5 ~ —30 MeV and B ~ 0.30, the experimental data on 
Q 2 F ril {Q 2 ) and Q 2 F rj i^{Q 2 ) can be consistently explained. It is found that fZ = —30 MeV 

Because we still have |/^;| << f Vc ~ 400 MeV, according to 



is consistent with Ref. 



41 
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FIG. 8: 7] - 7 (Left) and r/ — 7 (Right) transition form factors Q 2 F m (Q 2 ) and Q 2 F v > y (Q 2 ) with 



varying 1?, where 



-30 MeV, m n 



0.30 GeV, m s 



0.45 GeV and m r 



1.50 GeV. The 



experimental data are taken from Refs. 



13, 



33, 



36]. 



the mass-matrix-element shown by Ref. 



531 ] . we still have M Vc ~ m r r . up to high accuracy. 



And applying the parameters into the formulas presented in Ref. 53J, it can be found that 
the mixing between cc state with qq — ss basis is still quite small, i.e. the mixing is around 
1%. Then our present approximations (I17II18P are still reasonable. 



IV. SUMMARY 



Light pseudoscalar meson-photon transition form factor provides a good platform to study 
the leading-twist DA of the light pseudoscalar mesons since it contains only one bound state. 
In the present paper, we have analyzed three pseudoscalar meson-photon transition form 
factors consistently by using an uniform WF model suggested by Ref. 28j. By comparing 



the estimations with the experimental data on these form factors, it can provide strong 
constraints on the light pseudoscalar meson DAs. Our results are listed in the following. 

(1) According to Eqs.i PfTOlfTTl) . all pseudoscalar meson-photon transition form factors 
Q 2 F^(Q 2 ), Q 2 F riql {Q 2 ) and Q 2 F Vsl (Q 2 ) should have similar behaviors. Since no rapid 
growth of Q 2 F m (Q 2 ) and Q 2 F V ^(Q 2 ) as that of Q 2 F nJ (Q 2 ) has been found experimentally, 
these two form factors can be explained by setting B ~ 0.30 together with small charm 
quark component \ f^,\ ~ 30 MeV. Such a moderate DA with B ~ 0.30 for r] q and r] s , which 
corresponds to the second Gegenbauer DA moment around 0.35, may also be the pion DA 
behavior. 
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(2) We have made a detailed discussion on the form factors' uncertainties caused by 
the constituent quark masses m q and m s , the parameter B, the mixing angle and the 
possible intrinsic charm components f° and fZ. Firstly, the parameter B determines the 
main behavior of the form factors. By varying B G [0.00, 0.60], one can conveniently obtain 
the results for tt — 7 transition form factor with those obtained in the literature, in which 
DA behavior varies from AS-like to CZ-like accordingly. A CZ-like DA with B ~ 0.60 can 
explain the data in high Q 2 region, however it fails to explain the form factors' lower Q 2 
behavior. While the AS-like DA with B ~ 0.00 can give a better understanding for lower 
Q 2 region, it can not explain the present BABAR data for large Q 2 behavior. Secondly, 
the parameters m q , m s and <fi can improve the behavior slightly. With the increment of 
m q and m s , three form factors Q 2 F ni {Q 2 ) , Q 2 F rjl {Q 2 ) and Q 2 F rj i 1 {Q 2 ) increase in lower Q 2 
region and decrease in higher Q 2 region. Q 2 F m (Q 2 ) decreases with the increment of </>, 
while Q 2 F V > 1 (Q 2 ) increases with the increment of <fi. Thirdly, Q 2 F V1 (Q 2 ) and Q 2 F V / 1 (Q 2 ) 
are slightly affected by the charm component in low Q 2 region, while in high Q 2 region, the 
form factors are quite sensitive to fZ and they can be greatly suppressed by possible charm 
component. 

(3) It has been found that by adjusting these parameters within their reasonable regions, 
one can improve the estimations of the form factors to a certain degree but still can not 
solve the puzzle, especially to explain the behavior of 7r — 7 form factor within the whole 
Q 2 region consistently. Due to the cancelation between F Vql and F^, it is reasonable that 



Q 2 F r) i 1 {Q 2 ) tends to be a constant as Q 2 — > 00 [33j, |36]. However, it is hard to understand 
why Q 2 F m (Q 2 ) and Q 2 F 1T1 {Q 2 ) have such a quite different large Q 2 behavior. Especially 
it is hard to explain the rapid growth of Q 2 F n ^(Q 2 ), i.e. probably the logarithmic growth 
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22] , reported by BABAR collaboration |l4| to be consistent with the previously obtained 
low Q 2 behavior by CELLO and CLEO collaborations [12I. Possible charm components 
and fZ can shrink the gap between these two form factors to a certain degree, but it 
can not be the reason for such a large difference. We hope more experimental data on these 
form factors in large Q 2 region can clarify the present situation. If the BABAR data is 
confirmed, then there may indeed indicate new physics in these form factors, since it is hard 
to be explained by the current adopted light-cone pQCD framework. 
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